MEIXNER POLYNOMIALS IN SEVERAL VARIABLES 
SATISFYING BISPECTRAL DIFFERENCE EQUATIONS 
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Abstract. We construct a set 9Jtd whose points parametrize families of Meixner 
polynomials in d variables. There is a natural bispectral involution b on 931^ 
which corresponds to a symmetry between the variables and the degree in- 
dices of the polynomials. We define two sets of d commuting partial difference 
operators diagonalized by the polynomials. One of the sets consists of differ- 
ence operators acting on the variables of the polynomials and the other one on 
their degree indices, thus proving their bispectrality. The two sets of partial 
difference operators are naturally connected via the involution b. 

1. Introduction 

One of the reasons why the classical orthogonal polynomials p n (x) of a discrete 
variable x appear in numerous applications is the fact they are eigenfunctions of a 
second-order difference operator acting on x. A detailed account of their charac- 
teristic properties and application can be found in |16j . 

The last few years there have been advances in different aspects of the classifica- 
tion and the construction of orthogonal polynomials in several variables which are 
eigenfunctions of partial difference operators, see for instance [U H3 IH1 113 HU fT2j - 
The difficulty in the multivariate case is that the polynomials are no longer uniquely 
determined by the orthogonality measure (up to a multiplicative constant), and 
therefore their spectral properties depend on the way we apply the Gram-Schmidt 
process, see [3] for the general theory. 

Following the pioneering work of Duistermaat and Grunbaum [5] , we say that a 
family {P ni ,...,n 4 (%ii ■ ■ ■ > Xd)} of polynomials in d variables x±, . . . , Xd with degree 
indices tlx , . . . , rid solves a discrete-discrete bispectral problem, if there exist com- 
muting partial difference operators Cf,.. . acting on the variables xi, . . . ,x d 
with coefficients independent of n\ , . . . , rid an d commuting partial difference opera- 
tors £™, . . . , Cd acting on the degree indices n%, . . . , rid with coefficients independent 
of X\ , . . . , Xd such that 

£|P„ 1 ,...,„ (J (xi, . . . ,x d ) = Xj(ni, . . . ,n d )P ni ,...,n d (xi, . -.,x d ) (1.1a) 

£jP ni! ... <nd (xi, . ..,Xd) = Hj(xi, . . . , x d )P nu ... ind (xi,. . .,Xd), (1-lb) 

for j — 1, 2, . . . , d, where \j(n\, . . . ,rid) are independent ofxi, . . . ,x d and (J,j(x%, . . . , x d ) 
are independent of n±, . . . , n d . Within the context of orthogonal polynomials one 
is often interested in the commuting operators corresponding to multiplications by 
the independent variables Xj and therefore fj,j(xi, . . . , x d ) — Xj is one natural choice 
for the eigenvalues in (jl.lbj) . On the other hand, even in the one- variable case, the 
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eigenvalues in (|l.la|) can be quadratic functions (e.g. in the case of Hahn polynomi- 
als) and in higher dimensions, they can depend on some or all of the chosen degree 
indices n±,..., rid, see for instance the families discussed in (31 IS] ■ 

Recently |10j we proved that all families of multivariate Krawtchouk polynomi- 
als defined by Griffiths [5] are bispectral. We used a Lie-theoretic approach which 
naturally led to the construction of the partial difference operators in equations 
(|l.ip . Moreover, we showed how other known families of bispectral multivariate 
Krawtchouk polynomials fit within the above theory. The construction was a natu- 
ral extension of an earlier joint work with Terwilliger [11] for the bivariate polyno- 
mials defined in jSJ. We note that different proofs of the orthogonality properties 
of the multivariate Krawtchouk polynomials as well as interesting probabilistic ap- 
plications and connections to Gelfand pairs and character algebras were discussed 
in the papers Q31 EM HI] • 

In the present paper we give a similar parametrization of bispectral multivari- 
ate Meixner polynomials and the corresponding partial difference operators. The 
main result concerning the bispectrality of the multivariate Meixner polynomials 
constructed here is analogous to one in the Krawtchouk case [10] , but the proof is 
different. 

The paper is organized as follows. In the next section, we define a set 9Jt f j and for 
every m £ VJld and (3 £ C \ (—No) we construct multivariate Meixner polynomials 
{P ni ,...,n d (xi, ■ ■ ■ , xa; m, /3) : rij £ No}, mutually orthogonal on N$ with respect to 
a weight function of the form 



The construction is based on a generating function and is similar to one used by 
Griffiths [5] in the case of the Krawtchouk polynomials. In Section [3] we give 
a hypergeometric representation for the multivariate Meixner polynomials which 
yields the duality between the variables and the degree indices of the polynomials. 
The key ingredient is to "expand" the generating function appropriately and to 
regroup the terms by changing the summation indices. The computation is essen- 
tially the same as the one in the Krawtchouk case which can found in the work 
of Mizukawa and Tanaka [15], except that the expansion of the generating series 
now is infinite. In Section [4] we construct the commuting partial difference opera- 
tors {£j}j =1 , {£™}jLi and we prove that they are diagonalized by the polynomials 
{Pni,...,n d {x\, ■ ■ ■ ,Xd',vn., j5) : rij £ No} by using the generating function and the 
duality. Finally, we list two explicit families of Meixner polynomials in arbitrary 
dimension in Section [5] In the first example, we explain how we need to pick the 
point m £ 9Jld in order to obtain the multivariate Meixner polynomials defined 
in [12]. In the second example, we construct a point m £ 9Jld which depends on 
one free parameter. The corresponding Meixner polynomials can be thought of as 
analogs of the the multivariate Krawtchouk polynomials used in [1] . 




i=i 



2. Multivariate Meixner polynomials 



2.1. Basic notations and definitions. Throughout the paper we shall use stan- 
dard multi-index notations. For instance, if x = {x\, x%, . . . , Xd) £ C d and n = 
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(ni,ra2, • ■ • ,nd) £ Nq then 



•x^, n! = ni!n 2 ! ■ ■ -n d \, 



and |x| = xi + x 2 + ■ ■ ■ + Xd- 

First we introduce the set 9Jld whose points parametrize families of multivariate 
Meixner polynomials. In most applications, the parameters are real, but since 
the partial difference equations studied here have natural extensions for complex 
numbers, we work below with C rather than R. 

Definition 2.1. Let dK d denote the set of 4-tuples (cq,C,C,U), where cq is a 
nonzero complex number and C, C , U are (d + 1) x (d + 1) matrices with complex 
entries satisfying the following conditions: 

(i) C = diag(l, — ci, -c 2 , . . . , -c d ) and C* = diag(l, -ci, -c 2 . . . , — c^) are di- 
agonal; 

(ii) U — {ui,j)o<i,j<d is such that uqj — Uj t o = 1 for all j = 0, 1, . . . , d, i.e. 

(I 1 1 ... l\ 

I tii! U 12 ... 



17 = 



(2.1) 



(iii) The following matrix equation holds 

V t CVC = c Q I d+x , (2.2) 

where Z^+i denotes the identity {d + 1) X [d + 1) matrix. 

Note that Co is a nonzero number and therefore (12.21) implies that the matrices 
C, C and {/ in Definition 12.11 are invertible. We shall denote by c and c the d- 
dimensional vectors 

c = (ci,c 2 , . . . ,c d ), c = (ci,c 2 , . . . ,c d ). 

Comparing the entries on both sides of equation (|2.2[) for j = and j > 
we obtain the following identities: 

d 

^ c fc u fe! j = I - co5i,o (2.3a) 
fe=i 

2J CkUk,iUk,j — 1 + — <5ij. (2.3b) 
fc=i Cj 

Remark 2.2. We shall use later the involution b on dJld defined by 

fa : m = (c , C,C,U)->m= {co,C, C, U l ). (2.4) 

In particular, using (|2.3a[) with i — and fa we see that 

co = l-|c| = l-|c|. (2.5) 

Remark 2.3. If we start with nonzero numbers c\, C2, . . . , c<j such that cq = 1— |c| ^ 
we can define a point m = (c ,C,C,U) £ 971^ as follows. First, we construct a 
basis Vq,Vi, . . . ,Vd of column vectors in C d+1 mutually orthogonal with respect 
to the bilinear form (u>i,ui 2 ) = w\Cw2 such that vo — (1,1,.. .,1)*, the 0-th 
coordinate of Vj is 1 and (vj, Vj) ^ for j — 1,2, ... ,d- Then the matrix U with 
columns vq, v%, . . . , Vd is of the form given in (|2.ip and the diagonal matrix C is 
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uniquely determined by (|2.2|l . If d = 1, then v\ — (1, l/ci)* is uniquely determined 
from ci. However, when rf > 1 we have d(d — l)/2 degrees of freedom in choosing 
t/. 

In the rest of the paper we use {zj}^ =1 and {wj}j =1 to denote formal complex 
variables which are sufficiently close to 0. In all exponents below we fix the principal 
branch of the logarithmic function. 

For every m G DJld and /3 G C we consider the function of x — (x\, . . . ,Xd), 
z = (zi,...,Zd) 



G(x,z;m,0) = (l-\z\ 



-P—\ x \ 



nh-E 



(2.6) 



Definition 2.4. For m € 9Jt<j and f3 G C \ (—No) we define multivariate Meixner 
polynomials {P n (x; m, /?) : n G Nq} by expanding G(a;, z; m, /3) in a neighborhood 
of z = as follows 



G(x, z;m, /3) = V" ~-^- P n (x;m, 0)z 



(2.7) 



2.2. Orthogonality relations. Wc show next that the polynomials {P n (x; m, (3) : 
n G Nq} defined above are mutually orthogonal with respect to the weight (fi)\ x \c x / x\ 
on Nq. More precisely, using the notations introduced so far, the following theorem 
holds. 

Theorem 2.5. Suppose that \c\\ + |ca| + • • • + |q| < 1. Then for n,m G Nq we 

have 



^2 P n (x;m,P)P m (x;m,l3) 



(P)\x\ c x 



(2. 



Proof. Using the formula 



(1-1*1 



J2 (7) i fc i -fc 



/rt:<; 



fc! 



(2.9) 



wc find 



V G(x,z;m,/3)G(x,«;;m,/3)^ c I 
' a;! 



(l-HKi-M)-!^ i-£ 



-/J (2.10) 



If we set Z = (—1, zi, . . . , z^) 4 and IV = (— 1, Wi, . . . , Wd) 1 , then the expression in 
the big parentheses on the right-hand side above is equal to {UWyCiUZ). Using 
we find 



{UWfC{UZ) = c W t C- 1 Z = co 



ZjWj 



3 = 1 
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Thus, flUO} gives 

G(x,z;m,P)G{x,w;m,py-^ 



C = Cn 



ZjWj 



3 = 1 



The proof now follows by expanding the right side of the last equation using (12.9[) 
and comparing the coefficients of z n w m on both sides of the resulting identity. □ 

3. HYPERGEOMETRIC REPRESENTATION 

We denote by the set of all d x d matrices with entries in No . 
Theorem 3.1. We have 



A=(a ii3 )eM d 



(3.1) 

Note that (— nj)^ a . = when a « j n j- Therefore the sum in (|3.1[) is 

finite, using only matrices A with nonnegative integer entries satisfying a i,j — 

rij for j = 1, . . . , d. 

Proof of Theorem \3.1\ For every i e {1,2, we have 



3=1 



3 = 1 



Taking the Xi-th power of the above equation and expanding the second term using 
(12.91) we obtain 



'1,3^-3 



3 = 1 



= (i-Nr E 



(-l)l^l(-x t )| fci | 

fe! (l-|z|)l fc 'l 



where we have set ?7i = {1—u^x, 1 — 1^,2, 1 — Substituting the last formula 
into the right-hand side of (|2.6p we obtain 



G(x, z; m, /3) = 



J2 (i-^D-^-EtxI^I^^EtxI^I 
fei,/s2,...,fe<ieNg 

/ d 



n 

vi=l 



Expanding also (1 — \z\) 13 ^>=i l fci ' we get 

! Al ( ''+3:;-, i/- 

G(x, z; m, (3) 



l\ 



: E 

l,k x ,k 2 ,— ,k d e^ 



(3.2) 



n 
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Let us denote by m.j the entries of the vector ki, and by lj the entries of I, i.e. 
h = (a,i } i,ai,2, ■ ■ ■ ,a h d) and I = (h, h, ■ ■ ■ , Id)- Then 

/: " We?. - M (/, • v; ,„,.,): 

(3.3) 

The proof now follows by plugging (|3.3p in the first line of (|3.2j) and by replacing 
the sum over I with a sum over n = (n\, 712, . . . , Ud) where rtj = Zj + X)j=i a *j - '— ' 

As an immediate corollary of the hypergeometric representation (|3.ip and the 
involution b on 9Jtd in Remark 12.21 we see that the polynomials P n (x;m, fl) possess 
a duality between the variables {xj} and the indices {fij}. 

Corollary 3.2. For n,x £ Nq we have 

P n (x;m,p) = P x (n,m,l3). (3.4) 



4. BlSPECTRALITY 

Let {ei, e2, • ■ ■ , &d] be the standard basis for C d . For every i £ {1,2, ... ,d} 
we denote by E Xi and 2? n( the customary shift operators acting on functions of 
x = (x±, X2, ■ ■ ■ , a^d) and n = [n\, n%, ■ . ■ , nd), respectively, as follows 

E x J{x) = f(x + and E ni g(n) = g(n + e,). 

For fixed m S SDtd, /3 6 C and i 6 {1,2,..., d} we define the following difference 
operators 

= — Y"l CkUk,iUiiXi (E Xk E~ 1 - Id) 

d d 

- - V m,,^ - Id) - - V c fc u M (/3 + |x|) (£^ fe - Id) , (4.1a) 

and 

£? = — y c k u itk u ly ini (E nk E~J~ - Id) 

Cn * — ' 

l<kytl<d 

d d 

~ ~ u i,i n i (^"i 1 — Id) — — y CfcUi,fc(/3 + |n|)(_E„ - Id), (4.1b) 

c ° c ° is 

where Id denotes the identity operator. We show below that these operators are di- 
agonalized by the polynomials P n (x; m, /3), thus providing solutions to the bispectral 
problem. Note that the backward shift operators E~ x in (14.1b[) are multiplied by ni , 
and therefore these expressions will vanish when n; — 0. Thus Cf is a well-defined 
operator acting on functions of n S Nq. Similarly, in view of the orthogonality 
(|2.8|) . it is natural to consider the polynomials for x £ Nq and the operators Cf 
will involve evaluations of P n (a;;m, j3) only at points x with nonnegative integer 
coordinates. 
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Theorem 4.1. For every i G {1, 2, . . . , d} the following spectral equations hold 

C?P n (x; m, f3) = ni P n {x; m, (3), (4.2a) 
£?P n (x;m,p) =XiP n (x;m,P). (4.2b) 

Proof. Fix i G {1, 2, . . . , d} and let us denote 

= Id - B- 1 + ^ CfcU M (B^B- 1 - Id) for every Z = 1, 2, . . . , d 

fce{l,...,d}\{7} 

and 

d 

Ci 



Then 



Using the definition (|2 . 6[) of G(x, z; m, /?) we find 

G(x,z;m,P) £ W - Id) G^, z; m,^) 



fee{i,...,d}\{i} 
1 



(interchanging the sums) 



^ /XckUk,jUl,j ~ CkUk,iUk,j)Zj 



^ d d 



J 3 

(using equations (|2 - 3|) ) 

1 



1 - Ej=l j 3 = 1 fe=l 



^ /X c kUk,jUl,j - CkU k ,iU ktj )Zj 



i - E,=i «j 



3=1 V J 7 



J2(uij - l) Zj - j-Z t 



1 - Ei=i wjj^j; v=i 

From the last relation, the definition of and 



(Id - B x ) G(a;, z;m,/3) 



it follows that 

' £?,G(z,z;m,/3) = 7 ^ r . (4.4) 



A similar computation shows that 

g^-j^^^^ffl.^tWS. (4.5) 
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Using equations (|4.3|) . f)4.4[) and ()4.5p we see that 



£*G(x, z- m, 0) = z~G(x, z; m, /3), (4.6) 

which combined with (I2.7[) completes the proof of (|4.2ap . The proof of (|4.2b[) 
follows from the duality established in Corollary 13.21 Indeed, if a; G Nq, then 
P n (x;m, (i) = P x (n, m, /?) and equation (|4.2b[) follows from (|4.2a[) . Moreover, for 
fixed heNj both sides of (I4.2h>|) are polynomials in x of total degree at most \n\ + 1, 
hence if the equality holds for all x G Nq it will be true for arbitrary x G C d . □ 

5. Some examples 

In this section we illustrate how specific points in 9Jld lead to families of mul- 
tivariate Meixner polynomials or analogs of multivariate Krawtchouk polynomials 
which have already appeared in the literature in different applications. 

5.1. Let us fix nonzero complex numbers c\, . . . ,Cd and let Co — 1 — |c|. We define 
ci, C2, . . . , Cd in terms of ci, . . . , c<j as follows: 

5fe = ^ ~~~~ d f° r ^ = 1j 2, . . . , d. (5.1) 

(i -E J= fe+i c j)(i -Ej=fc c j) 

We shall assume that {cj} are such that the denominators in (|5.1[) do not vanish. 
Next we define a (d + 1) X (d+ 1) matrix U = (u^.j) with entries 

u i,j = <5o,i , when i < j, (5.2a) 

Mi.j = 1, when i > j, (5.2b) 

u ,o = 1, (5.2c) 

u t i = - 1 - ^fc^ 1 Cfc for i = 1,2,..., d. (5.2d) 

Thus, U is a matrix of the form (|2.1I) where the remaining entries are O's and 
l's above and below the diagonal respectively, and the diagonal entries are given 
in (|5.2d[) . One can check that with the above notations equation (|2 . 2[) holds and 
therefore we obtain a point m G 971^ which depends on the free parameters Ci, . . . , q. 
In particular, note that if G (0, 1) for all i G {0, 1, . . . , d} then Cj G (0, 1) for all 
i G {1, . . . , d}. Up to a permutation of the variables and the parameters, this choice 
leads to the multivariate Meixner polynomials defined in |12) . 

5.2. Fix now q ^ {0, 1} and define 

c fe = c k = (1 - q)q k - 1 for k = 1, 2, . . . , d. (5.3) 

Then Co = 1 — |c| — q d . We define also a (d + 1) x (d + 1) matrix U = {ui : j)o<i.j<d 
with entries 

Mi.j = 1, when i + j < d, (5.4a) 

Ui j — , when i + j = d + 1, (5.4b) 

1 — q 

Ui.j = 0, when i + j > d + 1. (5.4c) 

With the above notations, we see that (|2.2[) holds leading to a point m G 9Jt<j which 
depends on the free parameter q. Moreover, if q G (0, 1) then a = Ci G (0, 1) for all 
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i G {0, 1, . . • , d}. These Meixner polynomials can be thought of as analogs of the 
multivariate Krawtchouk polynomials used in [I] . 
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